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Using a custom 3 Cerenkov-ring fitter, we report cross sections for ^-induced charged-current 
single 7T° production on mineral oil (CH2) from a sample of 5810 candidate events with 57% signal 
purity over an energy range of 0.5 — 2.0 GeV. This includes measurements of the absolute total 
cross section as a function of neutrino energy, and flux-averaged differential cross sections measured 
in terms of Q 2 , fi~ kinematics, and 7r° kinematics. The sample yields a flux-averaged total cross 
section of (9.2 ± 0.3 s tat. ± 1.5 S3/s t.) X 10 -39 cm 2 /CH2 at mean neutrino energy of 0.965 GeV. 

PACS numbers: 13.15.+g, 25.30.Pt 



I. INTRODUCTION 

The charged-current interaction of a muon neutrino 
producing a single neutral pion (CC7r°) most commonly 



* Present address: California Institute of Technology; Pasadena, 

CA 91125 
t deceased 

t Present address: IFIC, Universidad de Valencia and CSIC; 46071 
Valencia, Spain 

§ Present address: Imperial College; London SW7 2AZ, United 
Kingdom 



occurs through the A(1232) resonance for neutrino ener- 
gies below 2 GeV. As there is no coherent contribution to 
CC7T° production, this process is an ideal probe of purely 
incoherent pion-production processes and thus offers ad- 
ditional kinematic information on 7r° production beyond 
what is measured in the neutral-current channel [l|, 0]. 
Previous measurements of CC7r° production at these en- 
ergies were made on deuterium at the ANL 12 ft bubble 
chamber Hi] and the BNL 7 ft bubble chamber @. To- 
tal cross-section measurements were reported on samples 
of 202.2 0] and 853.5 @ events for the ANL and BNL 
experiments respectively. Previous measurements 0] 
were also performed at higher neutrino energy on a vari- 
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ety of targets. 

Using Cerenkov light detection techniques, this pa- 
per revisits this topic and measures CC7r° production on 
carbon. In order to extract such interactions from the 
more dominant charged-current quasi-elastic (CCQE) 
and charged-current single 7r + (CC7T+) production pro- 
cesses, a custom fitter has been developed to isolate and 
fit both the fi~ and the 7r° in a CC71- event. This fitter 
also accurately reconstructs the kinematics of these in- 
teractions providing a means with which to extract both 
total and single-differential cross sections. Additionally 
reported is a measurement of the flux-averaged total cross 
section. This work presents the most comprehensive mea- 
surements of CC7T interactions to date, at energies be- 
low 2 GeV, on a sample of events 3.5 times that of the 
combined previous measurements. Results include the 
total cross section, the single-differential cross section in 
Q 2 , and the first measurements of single-differential cross 
sections in terms of final-state particle kinematics. The 
reported cross sections provide a combined measure of 
the primary interaction cross section, nuclear effects in 
carbon, and pion re-interactions in the target nucleus. 



II. FINAL STATE INTERACTIONS & 
OBSERVABLE CCtt 

Because this measurement is being performed on a nu- 
clear target, particular attention must be paid to how 
the sample is being defined, especially given how nuclear 
and final-state effects can influence the observables. The 
dominant effect is final state interactions (FSI) which are 
the re-interactions of particles created from the neutrino- 
nucleon interaction with the nuclear medium of the tar- 
get nucleus. FSI change the experimental signature of a 
neutrino-nucleon interaction. For example, if a n + from 
a CC7r + interaction is absorbed within the target nucleus 
and none of the outgoing nucleons are detected, then the 
event is indistinguishable from a CCQE interaction. Ad- 
ditionally, if the 7r + charge exchanges then the interaction 
is indistinguishable from a CC7r° interaction. This is due 
to the fact that the nuclear debris is typically unobserv- 
able in a Cerenkov-style detector. The understanding of 
FSI effects is model-dependent, with large uncertainties 
on the FSI cross sections. An "observable" interaction is 
therefore defined by the leptons and mesons that remain 
after FSI effects. Observable interactions are also inclu- 
sive of all nucleon final states. To reduce the FSI-model 
dependence of the measurements reported here, the sig- 
nal is defined as a \i~ and a single ir° that exits the target 
nucleus, with any number of nucleons, and with no ad- 
ditional mesons or leptons surviving the nucleus. This 
is referred to as an observable CC7r° event. The results 
presented here are not corrected for nuclear effects and 
intra-nuclear interactions. 



III. THE MINIBOONE EXPERIMENT 

The Mini Booster Neutrino Experiment (Mini- 
BooNE) Q is a high-statistics low-energy neutrino exper- 
iment located at Fermilab. A beam of 8 GeV kinetic en- 
ergy protons is taken from the Booster Q and impinged 
upon a 71 cm long beryllium target. The data set pre- 
sented in this paper corresponds to 6.27 x 10 20 p.o.t. (pro- 
tons on target) with an uncertainty of 2%. The result- 
ing pions and kaons are (de)focused according to their 
charge by a toroidal magnetic field created by an alu- 
minum magnetic focusing horn; positive charge selection 
is used for this data. These mesons then decay in a 50 m 
long air-filled pipe before the remnant beam impacts a 
steel beam dump. The predominantly f M -neutrino beam 
passes through 500 m of dirt before interacting in a spher- 
ical 800 ton, 12 m diameter, mineral oil (CH2), Cerenkov 
detector. The center of the detector is positioned 541 m 
from the beryllium target. The inner surface of the detec- 
tor is painted black and instrumented with 1280 inward- 
facing 8 inch photomultiplier-tubes (PMTs) providing 
11.3% photocathode coverage. A thin, optically-isolated 
shell surrounds the main tank region and acts to veto en- 
tering and exiting charged particles from the main tank. 
The veto region is painted white and instrumented with 
240 tangentially-facing 8 inch PMTs. A full description 
of the MiniBooNE detector can be found in Ref. [T(| . 

The neutrino beam is simulated within a Geant4 [ll| 
Monte Carlo (MC) framework. All relevant components 
of the primary proton beam line, beryllium target, alu- 
minum focusing horn, collimator, meson decay volume, 
beam dump, and surrounding earth are modeled [T3 | . 
The total p-Be and p-Al cross sections are set by the 
Glauber model Wherever possible, inelastic pro- 

duction cross sections are fit to external data. The neu- 
trino beam is dominated by produced by 7r + decay 
in flight. The 7T" 1 " production cross sections are set by 
a Sanford-Wa ng |14j ] fit to tt + production data provided 
by the HARP [15| and E910 [16[ experiments. The high- 
energy {E v > 2.4 GeV) neutrino flux is dominated by 
from K + decays. The K + production cross sections 
are set by fitting data from Refs. [l7]-[24j to a Feynman 
scaling parametrization. The production of protons and 
neutrons on the target are set using the MARS [25| sim- 
ulation. The v e , and v e contributions to the flux are 
unimportant for this measurement. Ref. [l2j describes 
the full details of the neutrino-flux prediction and esti- 
mation of its systematic uncertainty. It should be noted 
that the MiniBooNE neutrino data has not been used to 
tune the flux prediction. 

Interactions of neutrinos with the detector materials 
are simulated using the v3 Nuance event generator [26( . 
The Nuance event generator is a comprehensive simu- 
lation of 99 neutrino and anti-neutrino interactions on 
nuclear targets over an energy range from 100 MeV 
to 1 TeV. The dominant interaction in MiniBooNE, 
CCQE, is modeled according to Smith- Moniz [27J; how- 
ever, the axial mass, Ma, has been adjusted for better 
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agreement with the MiniBooNE data to M'f = 1.23 ± 
0.20 GeV/c 2 (H. The target nucleus is simulated with 
nucleons bound in a relativistic Fermi gas (RFG) [13] 
with binding energy Eb — 34 ± 9 MeV and Fermi mo- 
mentum p F = 220 ± 30 MeV/c ^ (on carbon). The 
RFG model is further modified by shape fits to Q 2 , for 
bette r ag reement of the CCQE interaction to MiniBooNE 
data [Hf . The Rein-Sehgal model [3(| is used to predict 
the production of single-pion final states for both CC and 
NC modes. This model includes 18 non-strange baryon 
resonances below 2 GeV in mass and their interference 
terms. The model [3(| predicts the A(1232) resonance 
to account for 71% of CC7r° production (84% of resonant 
production); a 14% contribution from higher mass reso- 
nances; a 15% contribution from non-resonant processes. 
The non-resonant processes are added ad hoc to improve 
agreement to past data 30]. They and are not indica- 
tive of actual non-resonant contributions but allow for 
additional inelastic contributions. The quarks are mod- 
eled as relativistic harmonic oscillators a la the Feynman- 
Kislinger-Ravndal model [3l[ . The axial mass for single- 
pion production is Mjf = 1.10 ± 0.27 GeV/c 2 , and for 
multi-pion production it is = 1.30 ± 0.52 GeV/c 2 . 

The model includes the re-interactions of both baryon 
resonances, pions, and nucleons with the spectator nu- 
cleons leading to the production of additional pions, pion 
charge-exchange, and pion absorption. For observable 
CC7T+ interactions, the model is reweighted to match the 
MiniBooNE data by a technique described in §VI Al 

The MiniBooNE detector is simulated within a 
Geant3 [32| MC. This MC handles the propagation of 
particles after they exit the neutrino-target nucleus, sub- 
sequent interactions with the mineral oil [331 ]. and most 
importantly, the propagation and interactions of optical 
photons. Photons with wavelengths between 250-650 nm 
are considered. The production, scattering, fluorescence, 
absorption, and reflections of these photons are mod- 
eled in a 35 parameter custom optical model [l^, H3, HI] . 
The detector MC also simulates photon detection by the 
PMTs and the effects of detector electronics. The ab- 
sorption (tt + \) and charge- exchange (ir + — > tt°) of 
7r + particles on carbon are fixed to external data (36l - [38j 
with uncertainties of 35% and 50% respectively. 



IV. EVENT RECONSTRUCTION 

Particles traversing the mineral oil are detected by the 
Cerenkov and scintillation light they produce. The rela- 
tive abundances of these emissions, along with the shape 
of the total Gerenkov angular distribution, are used to 
classify the type of particle in the detector. For a single 
particle, an "extended track" is fit using a maximum- 
likelihood method for several possible particle hypothe- 
ses. For each considered particle type, the likelihood is 
a function of the initial vertex, kinetic energy, and direc- 
tion. For a given set of track parameters, the likelihood 
function calculates probability density functions (PDFs) 



for each of the 1280 PMTs in the main portion of the de- 
tector [H, HI] . Separate PDFs for the initial hit time and 
total integrated charge are produced. As the data acqui- 
sition records only the initial hit time and total charge 
for each PMT, the Gerenkov and scintillation contribu- 
tions are indistinguishable for a given hit; however, they 
are distinguishable statistically. The likelihood is formed 
as the product of the probabilities, calculated from the 
PDFs. The initial track parameters are varied using Ml- 
nuit [40j and the results of the best fit likelihood de- 
termines the parameters for both the particle type and 
kinematics. 

The extended-track reconstruction is scalable to any 
number of tracks. The charge PDFs are constructed by 
adding the predicted charges from each track together to 
determine the overall predicted charge. The time PDFs 
are calculated separately for each track, and separately 
for the Cerenkov and scintillation portions, then time 
sorted and weighted by the probability that a particu- 
lar PDF caused the initial hit [TH [34| . The reconstruc- 
tion needed for an observable CC7r° event requires three 
tracks: a ix~ track and two photon tracks from a common 
vertex [4lj, [42|. The final state is defined by a a tt°, 
and nuclear debris. The [i~ is directly fit by the recon- 
struction, along with its decay electron. The 7r° decays 
into two photons with a branching fraction of 98.8% [43| 
at the neutrino interaction vertex (ct = 25.1 nm). The 
two photons are fit by the reconstruction. Photons pro- 
duce Cerenkov rings both by converting (A = 67 cm) into 
e + e~ pairs through interactions with the mineral oil and 
by Compton scattering. The nuclear debris is ignored in 
the reconstruction as it is rarely above Cerenkov thresh- 
old and therefore only contributes to scintillation light. 
As the calculation of the kinetic energy of a track is dom- 
inated by the Cerenkov ring, the added scintillation light 
is effectively split uniformly among the three tracks and 
justifiably ignored. 

The novelty of the CC7r° event reconstruction is the 
ability to find and reconstruct three Cerenkov rings. The 
CC7T likelihood function is parametrized by the event 
vertex (a;, y, z), the event start time (t), the /i~ direction 
and kinetic energy (d^, (p^, E^), the first 7 direction and 
energy (di,<px,E\), the first 7 conversion length (si), the 
second 7 direction and energy (#2, 4>2, E2), and conver- 
sion length (S2). Ring, or track finding is performed in 
a stepwise fashion. The first track is seeded in the likeli- 
hood function and fixed by the /j, one-track fit described 
in Ref. [39[. The second track is scanned through 400 
evenly-spaced points in solid angle assuming 200 MeV 
kinetic energy and no 7 conversion length. The best scan 
point is then allowed to float for both tracks simultane- 
ously. The third track is found by fixing the two tracks 
and scanning again in solid angle for the third track. 
Once the third track is found, the best-scan point, along 
with the two fixed tracks, are allowed to float. This stage 
of the fit, referred to as the "generic" three- track fit, de- 
termines a seed for the event vertex, track directions, 
and track energies. A series of three parallel fits are then 



4 



performed to determine track particle types and the fi- 
nal fit kinematics. Each of these fits is seeded with the 
generic three-track fit. The fits assign a /i~ hypothesis 
to one of the tracks and a 7 hypothesis to the other two, 
allowing for the possibility that the [i~ was not found 
by the original /1 one-track fit and was found during the 
first or second scan. The conversion lengths are seeded at 
50 cm and fit along with the kinetic energies while keep- 
ing all the other parameters fixed, thereby determining 
their seeds for the final portion of the fit. The final stage 
of each fit allows all parameters to float, taking advantage 
of Minuit's Improve function [iol |. A term is added to 
the fit negative-log-likelihoods comparing the direction 
from the fit event vertex to the fit pT decay vertex ver- 
sus the fit /i~ direction weighting by the separation of the 
vertices. This term improves the identification of the par- 
ticle types. The likelihoods are then compared to choose 
the best fit. For further details see Ref. [41j . 
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FIG. 1. Top: Fractional Gaussian resolutions as functions 
of T M (solid) and |p^o|c (dashed). Bottom: Angular resolu- 
tions as functions of cos 6^ (solid) and cos^o (dashed). The 
resolutions between each particle type are also correlated. 

The quality of the reconstruction is assessed by evalu- 
ating the residual resolutions of the signal using the MC. 
Fig. Q] shows the residual resolutions for both the p~ and 
the 7T . The overall fi~ kinetic-energy fractional resolu- 
tion is 7.4%, and the angular resolution is 2°. The 7r°, 
being a combination of the two fit photons, has an overall 
momentum resolution of 12.5% and angular resolution of 
7.8°. The energy and momentum resolutions are worse 
at lower energies and momenta and flatten out toward 
larger values. The \i~ angular resolution is mostly flat 
and is only slightly better in the forward direction. The 
7r° angular resolution gets much better in the forward 
direction as forward going ir° tend to have larger mo- 



mentum. The resolutions between each particle type are 
also somewhat correlated. Additionally, the interaction- 
vertex resolution is 16 cm. 

The initial neutrino energy is calculated, assuming that 
the signal events are from the reaction v^n — > fi~ir°p, 
from the measured yT and ir° kinematics under three as- 
sumptions: the interaction target is a stationary neutron, 
the hadronic recoil is a proton, and the neutrino is travel- 
ing in the beam direction. Under these assumptions, E v 
is constrained even if the proton is unmeasured [42| . The 
assumption of a stationary neutron contributes to smear- 
ing of the reconstructed neutrino energy because of the 
neutron's momentum distribution. The neutrino-energy 
resolution is 11%. Ideally one would measure the pro- 
ton, and additional hadronic debris; however, as these 
particles are rarely above Cerenkov threshold, it is im- 
practical to do so in MiniBooNE. Also, as only 70% of 
the observable CC7r° interactions are nucleon-level CC7r° 
on neutrons, additional smearing is due to CC7r + charge- 
exchanges on protons or other inelastic processes produc- 
ing a 7T° in the final state. These smearings are not ex- 
pected to be large. The 4-momentum transfer, Q 2 , to the 
hadronic system can be calculated from the reconstructed 
muon and neutrino momentum (42j. The calculation of 
the nucleon resonance mass is performed from the neu- 
trino and muon momentum also assuming a stationary 
neutron (see appendix). 



V. EVENT SELECTION 

Isolating observable CC7r° interactions is challenging 
as such events are expected to comprise only 4% of the 
data set [26[ . The sample is dominated by observable v^- 
CCQE (44%), with contributions from observable CCn + 
(19%), and other CC and NC modes. Basic sorting is first 
performed to separate different classes of events based on 
their PMT hit distributions. Then the sample is further 
refined by cutting on reconstructed quantities to yield an 
observable CC7r° dominated sample. Each cut is applied 
and optimized in succession and will be discuss over the 
remainder of this section. 

The detector is triggered by a signal from the Booster 
accelerator indicating a proton-beam pulse in the Booster 
neutrino beam line. All activity in the detector is 
recorded for 19.2 /js starting 4.6 fj,s before the 1.6 fj,s 
neutrino-beam time window. The detector activity is 
grouped into "subevents:" clusters in time of PMT hits. 
Groups of 10 or more PMT hits, within a 200 ns win- 
dow, with time spacings between the hits of no more 
than 10 ns with at most two spacings less than 20 ns, 
define a subevent jic| . In an ideal neutrino event, the 
first subevent is always caused by the prompt neutrino 
interaction; subsequent subevents are due to electrons 
from stopped-muon decays. Neutrino events with one 
subevent are primarily due to neutral-current interac- 
tions and j/ e -CCQE. Two-subevent events are from v^- 
CCQE, CC7T , and NC7r + . Three-subevent events are al- 
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most completely CC7r + , with some multi-7r production. 
Stopped muon decays produce electrons with a maximum 
energy of 53 MeV which never cause more than 200 PMT 
hits in the central detector. 

To select a sample of contained events, a requirement 
is made to reject events that penetrate the veto. These 
events typically cause more than 6 veto PMT hits. There- 
fore, a two-subevent sample is defined by requiring more 
than 200 tank PMT hits in the first subevent, fewer than 
200 tank PMT hits in the second, and fewer than 6 veto 
PMT hits for each subevent. The two-subevent sam- 
ple is predicted to be 71% z^-CCQE, 16% CCtt + , and 
6% CC7T . Observable CC7r + events make it into the 
two-subevent sample for several reasons: primarily by 
7r + — > 7T° and 7r + — > ^ in mineral oil, by fast muon de- 
cays whose electrons occur during the prior subevent, and 
by /i~ capture on nucleons affecting 8% of /i~ in mineral 
oil. The two-subevent filter keeps 40% of CCQE and 
CC7T interactions while rejecting 80% of CC7r + interac- 
tions. The rejection of CCQE and signal is mainly from 
/i~ that exit the tank. Additionally, events are lost by 
the 8% /i~ capture rate on carbon. 
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FIG. 2. (color online) The one-track fit likelihood ratio vs. the 
one-track muon fit kinetic energy. A separation of CCQE 
events (red "o") relative to CC7r° events (grey "x") is per- 
formed by selecting events below the black line. For clarity, 
the events plotted are prescaled by 1000 and 100 for CCQE 
and CC7T respectively. The cut is optimized on the full non- 
prescaled MC sample. 

To isolate a purer sample of observable CC7r° events 
from the two-subevent sample, before the observable 
CC7T fit is performed, ^ M -CCQE events are rejected by 
cutting on the ratio of the f^-CCQE fit likelihood to the 
j/ e -CCQE fit likelihood as a function of //"-fit kinetic en- 
ergy (see Fig. [5]). This cut is motivated by the fact that 
i^-CCQE interactions are dominated by a sharp muon 
ring, while CC7T interactions (with the addition of two 



photon rings) will look "fuzzier" and more electron-like 
to the fitter. The cut, which has been optimized to re- 
ject CCQE, rejects 96% of observable j^-CCQE while 
retaining 85% of observable CC7r° [4l[. Additionally, a 
reconstructed radius cut of r rcc < 550 cm is used to con- 
strain events to within the fiducial volume. 
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FIG. 3. The smallest angle between two of the three recon- 
structed tracks. Displayed are the data (with statistical er- 
rors), total MC (solid), identified fi~ (dashed), and mis-IDed 
fi~ (dotted). The cut selects events above 0.6 radians. 

To reject misreconstructed signal events, along with 
certain backgrounds, a cut is applied to reject events if 
two of the three Cerenkov rings reconstruct on top of 
one another. When this occurs the fitter ambiguously 
divides the total energy between the two tracks. For 
cases where there are two or fewer rings, the fitter will 
place two of the tracks in the same direction. For cases 
where there are three or more Cerenkov rings, the fit- 
ter can still get trapped with two tracks on top of one 
another. This can happen either because of asymmet- 
ric 7T° decays, a \i~ near or below Cerenkov threshold, a 
dominant fi~ ring, or events that truly have overlapping 
tracks. In all of these cases, the reconstruction becomes 
poor, especially the identification of the fi~ in the event. 
Fig. [3] shows the smallest reconstructed angle between 
two of the three reconstructed tracks for both data and 
MC. The MC events are separated into samples that cor- 
rectly identified the /i~ and those that did not. A cut is 
optimized on signal events in the MC to reject misiden- 
tified fi~ , and rejects events with track separations less 
than 0.6 radians [4l|. This cut reduces the expected \x~ 
misidentification rate to the 20% level. Additionally the 
observable CC7r° fraction is increased to 38% and the 
observable CCQE fraction is reduced to 13%. 

The next series of cuts reject non-7r° backgrounds. The 
first requirement compares the observable CC7T fit like- 
lihood vs. a generic three-track fit and selects events that 
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FIG. 4. The logarithm of the ratio of the observable CC7r° fit 
likelihood to a generic three-track fit. Displayed are the data 
(with statistical errors), total MC (solid), observable CC7r° 
(dashed), backgrounds with a 7r° in the final state or produced 
after the event (dotted), and backgrounds with no n° (dot- 
dashed). The cut selects events above 0.06. 



FIG. 5. The reconstructed 77 mass. Displayed are the data 
(with statistical errors), total MC (solid), observable CC7r° 
(dashed), backgrounds with a TV in the final state or produced 
after the event (dotted), and backgrounds with no n° (dot- 
dashed). The vertical dotted line is the known tv° mass. 



are more CC7r°-like. As the observable CC7r° fit is ac- 
tually a generic-^77 from a common vertex-fit, this cut 
selects events that match this criterion. The second cut 
is on the reconstructed 77 mass about the expected 7T° 
mass. This cut demands that the photons are consis- 
tent with a 7T° decay. The combination of these cuts 
define the observable CC7r° sample. Fig.S]shows the log- 
arithm of the ratio of the observable CC7T fit over the 
generic three-track fit likelihoods. The MC is separated 
into three samples: observable CC7r°, background events 
with a 7T° in the final state or created later in the event, 
and background events with no ir°. Both the observ- 
able CC7T and backgrounds with a ir° are more /ijj-like 
than events with no ir° in the event. Additionally, as the 
backgrounds with a tt° either have multiple pions or the 
7r was produced away from the event vertex, the likeli- 
hood ratio for these events tend slightly more toward the 
generic fit. Events with no ir° peak sharply at low ratio 
values. The optimization rejects non-7r° backgrounds by 
selecting events greater than 0.06 in this ratio 41 1. 

The final cut on the reconstructed 77 mass defines the 
observable CC7r° event sample. Fig. [5] shows the recon- 
structed 77 mass for both data and MC. No assump- 
tion is used in the fit that the two photons result from 
a 7T° decay; nevertheless, both data and MC peak at the 
known ir° mass. The predicted background MC with a 
7T° in the final state, or a 7r° produced after the event, 
has a broader peak than the signal MC. This broadening 
occurs for the same reasons discussed for the likelihood 
ratio; these events either produced a tt° away from the 
/i - vertex, or there are multiple pions in the final state. 



TABLE I. The expected efficiency and purity of observable 
CC7r° events as a function of applied cut. 



cut description 



efficiency purity 



none 100% 3.6% 

Two-subevent and Tank and Veto hits 38.2% 5.6% 

CCtt filter and fiducial volume 27.9% 29.6% 

Misreconstruction 10.3% 38.1% 

Likelihood ratio and 6.4% 57.0% 



As one might expect, background events with no ir° in 
the final state show no discernible mass peak and pile up 
at low mass with a long misreconstructed tail extending 
out to high mass. A cut is optimized to select events 
around the known tt° mass [0.09 < m 77 < 0.2 GeV] to 
reject non-7r° backgrounds (low mass cut) and to increase 
signal purity (high mass cut) [4l| . The addition of these 
selection cuts increases the observable CC7r° purity to 
57% with 6% efficiency. After all cuts, the observable 
CC7T candidate sample contains 5810 events in data for 
6.27 x 10 20 p.o.t. while the MC predicts 4160.2 events. 
Table Q] summarizes the effects of the cuts on the MC 
sample, while Table |TT] summarizes the background con- 
tent of the observable CC7r° candidate sample. 



VI. ANALYSIS 

The extraction of observable CC7r° cross sections from 
the event sample requires a subtraction of background 
events, corrections for detector effects and cut efficiencies, 
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TABLE II. The expected background composition of the 
CC7T candidate sample by observable mode. The level of 
observable CC7r + events are determined by the method de- 
scribed in Section TVI Al The symbol X represents all nuclear 
final states and photons, and N > 2 and M > 1 are the num- 
ber of pions in the final state. Other backgrounds include 
deep inelastic scattering and NC elastic scattering. 



observable 






fraction of 


mode 


description 


background 


CC7T+ 


«/„CH 2 - 


-> fl-TT+X 


52.0% 


CCQE 


VpCHa - 


->■ pTX 


15.4% 


CCmulti-7r 


f M CHa - 


-> ^-(Nti)X 


14.0% 


NC-TT 


v CH 2 - 


■¥ v (Mtv)X 


8.8% 


others 






9.8% 



a well-understood flux, and an estimation of the number 
of interaction targets. The cross sections are determined 
by 



da 

dx 



n<Fe, Air, 



(1) 



where x is the variable of interest, i labels a bin of the 
measurement, Axi is the bin width, Nj is the number of 
events in data of bin j, Bj is the expected background, 
Uij is a matrix element that unfolds out detector effects, 
6j is the bin efficiency, $j is the predicted neutrino flux, 
and n is the number of interaction targets. For the single- 
differential cross-section measurements, the flux factor, 
$i, is constant and equals the total flux. For the total 
cross-section measurement as a function of neutrino en- 
ergy the flux factor is per energy bin. The extracted cross 
sections require a detailed understanding of the measure- 
ments or predictions of the quantities in Eqn. [1] and 
their associated systematic uncertainties. By construc- 
tion, the signal cross-section prediction has minimal im- 
pact on the measurements. Wherever the prediction can 
affect a measurement, usually through the systematic er- 
ror calculations, the dependence is duly noted. 



A. CC7r + backgrounds 



ground, and tt + — >• ir° and ir + — > \ processes in the min- 
eral oil have large uncertainties. In particular, ir + 7r° 
in mineral oil is, by definition, not an observable CC7r° 
because the tt° did not originate in the target nucleus. 
By tying the observable CCn + production to measure- 
ments within the MiniBooNE data, the uncertainty on 
this background can also be further reduced. The to- 
tal error is separated into an uncertainty on CC7r + pro- 



duction and an uncertainty on tt 



and tt 



X 



processes in mineral oil occurring external to the initial 
target nucleus. Using the high statistics MiniBooNE 3- 
subevent sample, many measurements of the absolute ob- 
servable CCtt + production cross sections have been per- 
formed [iij Hfij]. This observable CC7r + sample is pre- 
dicted to be 90% pure [III |4(| making it the purest mode 
measured in the MiniBooNE data set to date. It is a 
very useful sample for this analysis as the bulk of the 
CC-7r + background events in the CC7r° sample fall in the 
kinematic region well-measured by the observable CCn + 
data. The measurements used are from the tables in 
Ref. [46|. Because the re- interaction cross sections for 
tt + — > 7T° and 7r + — > )r are strong functions of ir + energy, 
the constraint on CC7T + events is applied as a function of 
7r + kinetic energy and neutrino energy. Fig. [6] shows the 




The first stage of the cross-section measurement is to 
subtract the expected background contributions from the 
measured event rate. This is complicated by the fact that 
previously measured modes in MiniBooNE (CCQE (44j |. 
NC7T [![, and CC7T + (45|) show substantial normaliza- 
tion discrepancies with the Nuance prediction. The sin- 
gle largest background, observable CC7r + , is well con- 
strained by measurements within the MiniBooNE data 
set. The observable CCQE in the sample is at a small 
enough level not warrant further constraint. Most of the 
remaining backgrounds are unmeasured but individually 
small. 

The CC7r + backgrounds are important to constrain 
for two reasons: they contribute the largest single back- 



FIG. 6. (color online) The ratio of the measured observable 
CC7r + cross section to the Nuance prediction. The ratio is 
plotted as a function of and E v . 

ratio of the measured observable CC7T+ cross section as a 
function of 7r + and neutrino energies to the Nuance pre- 
dicted cross section. As the differential cross section with 
respect to ir + kinetic energy was not reported for bins 
with small numbers of events, the reweighting function is 
patched by the ratio of the total observable CC7T+ cross 
section as a function of neutrino energy. This reweighting 
factor (as high as 1.6 in some bins) is applied to every 
observable CC7r + event in the MC sample. All figures 
and numbers presented in this article have this reweight- 
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ing applied. By using the MiniBooNE data to constrain 
the CC7r + backgrounds, strict reliance on the Nuance 
implementation of the Rein-Sehgal model to predict this 
important background is avoided. 



B. Background subtraction 

The MC predicts a sample that is 57% pure observ- 
able CC7T after all analysis cuts. A 23% contribution 
of observable CC7r + interactions is set by prior measure- 
ments in the MiniBooNE data. The remaining 20% of the 
event sample is mostly comprised of CCQE and multi- 
7r final-state events (see Table Hi]) . These backgrounds, 
while thought to be produced in larger quantities than 
the MC prediction, are set by the MC as there is no clear 
method to extract the normalization of many of these 
modes from the current data. If the normalizations are 
diffcrcnt-but on the same order as previously measured 
modes-they would change the final results by at most a 
few percent. The uncertainties applied to the produc- 
tion of these backgrounds more than cover the possible 
normalization differences. After subtracting the back- 
grounds from the data the sample contains 3725.5 signal 
events with E v 6 0.5 — 2.0 GeV, which should be com- 
pared to 2372.2 predicted events. This is a normalization 
difference of 1.6, the largest normalization difference that 
has been observed in the MiniBooNE data thus far. This 
is attributed to the large effects of pion re-interactions 
which can directly influence the measured cross section 
for observable CC7r° production and the fact that the 
Nuance prediction appears low even when compared to 
prior measurements on deuterium data (47| . 



C. Unfolding & flux restriction 

To correct for the effects of detector resolution and 
reconstruction, a method of data unfolding is per- 
formed [48] . The unfolding method constructs a response 
matrix from the MC that maps reconstructed quantities 
to their predicted values. The chosen method utilizes the 
Bayesian technique described in Ref. 49]. This method 
of data unfolding requires a Bayesian prior for the sig- 
nal sample, which produces an intrinsic bias. This bias 
is the only way that the signal cross-section model af- 
fects the measured cross sections. In effect, this allows 
the signal cross-section model to pull the measured cross 
sections toward the shape of the distributions produced 
from the model while preserving the normalization mea- 
sured in the data; nevertheless, strict dependence on the 
signal cross-section is avoided. For these measurements, 
the level of uncertainty on the signal model cover the ef- 
fect of the bias on the central value. For situations where 
the unfolding is applied to a distribution that is signifi- 
cantly different than the Bayesian prior, the granularity 
of the unfolding matrix needs to be increased to stabi- 
lize the calculation of the systematic uncertainty. Oth- 



erwise the uncertainties would be larger than expected 
due to larger intrinsic bias. For the total cross sec- 
tion as a function of neutrino energy, the unfolding acts 
on the background-subtracted reconstructed neutrino en- 
ergy and unfolds back to the neutrino energy prior to 
the interaction. For each of the flux-averaged differen- 
tial cross-section measurements, the unfolding acts on 
the two-dimensional space of neutrino energy and the re- 
constructed quantity of a particular measurement. For 
final-state particles, /Lt~ and 7r°, the unfolding corrects to 
the kinematics after final-state effects, and are the least 
model dependent. For the 4-momentum transfer, Q 2 , 
the unfolding extracts to Q 2 calculated from the initial 
neutrino and the final-state \i~ . In both the total cross 
section and flux-averaged differential cross section in Q 2 , 
the final-state interaction model does bias the measure- 
ments. Uncertainties in the signal cross-section model, 
along with the final-state interaction model, arc expected 
to cover this bias. The unfolded two-dimensional distri- 
butions are restricted to a region of unfolded neutrino 
energy between 0.5 — 2.0 GeV. This effectively restricts 
the differential cross-section measurements over the same 
range of neutrino energy, and flux, as the total cross- 
section measurement. 



D. Efficiency correction 

The unfolded distributions are next corrected by a bin- 
by-bin efficiency. The efficiencies are estimated by taking 
the ratio of MC signal events after cuts to the predicted 
distribution of signal events without cuts but restricted 
to the fiducial volume. The efficiency is insensitive to 
changes in the underlying MC prediction to within the 
MC statistical error. While this additional statistical er- 
ror is not large, it is properly accounted for in the error 
calculation. The overall efficiency for selecting observ- 
able CC7T interactions is 6.4%; the bulk of the events 
(~ 60%) are lost by demanding a \i~ that stops and de- 
cays in the MiniBooNE tank. The cuts to reduce the 
backgrounds and preserve well-reconstructed events ac- 
count for the remainder. 



E. Neutrino flux 

A major difficulty in extracting absolute cross sections 
is the need for an accurate flux prediction. The Mini- 
BooNE flux prediction [l2j comes strictly from fits to ex- 
ternal data fla-KMj and makes no use of the MiniBooNE 
neutrino data. Fig. shows the predicted flux with 
systematic uncertainties. The flux is restricted by the 
unfolding method to the range 0.5-2.0 GeV, although 
contributions from the flux outside of this range affect 
the systematic uncertainties of the unfolding method. 
The integrated flux over this range is predicted to be 
= (3.545 ±0.259) x lO" 10 j/ M /p.o.t./cm 2 . 
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FIG. 7. (color online) The predicted flux with systematic 
uncertainties over the range 0.5-2.0 GeV. A table containing 
the values of the flux, along with a full-correlation matrix, is 
provided in Table. IVl of the appendix. 



F. Number of targets 

The Marcol 7 mineral oil is composed of long chains 
of hydrocarbons. The interaction target is one link in a 
hydrocarbon chain; CH2 chains with an additional hy- 
drogen atom at each end. The molecular weight is the 
weight of one unit of the chain averaged over the av- 
erage chain length. The density of the mineral oil is 
Ami = 0.845 ±0.001 g/cm 3 [io| . Thermal variations over 
the course of the run were less than 1%. The fiducial vol- 
ume is defined to be a sphere 550 cm in radius. Therefore, 
there are (2.517 ± 0.003) x 10 31 interaction targets. 



VII. SYSTEMATIC UNCERTAINTIES 

Sources of systematic uncertainties are separated into 
two types: flux sources and detector sources. Flux 
sources affect the number, type, and momentum of the 
neutrino beam at the detector. Detector sources affect 
the interaction of neutrinos in the mineral oil and sur- 
rounding materials, the interactions of the produced par- 
ticles with the medium, the creation and propagation of 
optical photons, and the uncertainties associated with 
the detector electronics. Whenever possible, the input 
sources to the systematic uncertainty calculations are fit 
to both ei-and in-situ data and varied within their full 
error matrices. These variations are propagated through 
the entire analysis chain and determine the systematic er- 
ror separately for each quantity of interest. For sources 
that affect the number of interactions (e.g. flux and cross- 
section sources), the central value MC is reweighted by 



taking the ratio of the value of the underlying parame- 
ter's excursion to that of the central value. This is done 
to reduce statistical variations in the underlying system- 
atics. Sources that change the properties of an event 
require separate sets of MC to evaluate the error matri- 
ces. The unfortunate aspect of this method is it intro- 
duces additional statistical error from the excursion. If 
the generated sets are large, then this additional error is 
small. 



A. Flux sources 

The main sources of flux uncertainties come from parti- 
cle production and propagation in the Booster neutrino 
beam line. The proton-beryllium interactions produce 
7r + , 7r~, K + , and K° particles that decay into neu- 
trinos. The dominant source of in the flux range of 
0.5-2.0 GeV is n + decay. The n + production uncertain- 
ties are determined by propagating the HARP measure- 
ment error matrix by spline interpolation and reweight- 
ing into an error matrix describing the contribution of 7r + 
production to the MiniBooNE flux uncertainty (l2^ . 
Over the flux range of this analysis, the total uncertainty 
on 7r + production introduces a 6.6% uncertainty on the 
i/^-flux. These uncertainties would be larger if the en- 
tire neutrino flux were considered in these measurements; 
however, as the largest flux uncertainties are in regions 
of phase space that cannot produce a CC7r° interaction, 
it was prudent to restrict the flux range. Horn related 
uncertainties stem from horn-current variations and skin- 
depth effects (which mainly affect the high-energy neu- 
trino flux), along with other beam related effects (e.g. 
secondary interactions), and provide an additional 3.8% 
uncertainty. No other source contributes more than 0.2% 
(K + production) to the uncertainty. The total uncer- 
tainty on the flux over the range 0.5-2.0 GeV, from all 
flux-related sources, is 7.3%. As the flux prediction af- 
fects the measured cross sections through the flux weight- 
ing, background subtraction, and unfolding, the result- 
ing uncertainties on the measured total cross section are 
7.5% and 7.3% for horn variations and tt + production 
respectively. 



B. Detector sources 

Uncertainties associated with the detector result from: 
neutrino-interaction cross sections, charged-particle in- 
teractions in the mineral oil, the creation and interaction 
of photons in the mineral oil, and the detector read-out 
electronics. The neutrino-interaction cross sections are 
varied in Nuance within their error matrices. These 
variations mainly affect the background predictions; how- 
ever, the expected signal variations (i.e. flux and signal 
cross-section variations) do affect the unfolding. The 
variations of the background cross sections cover the dif- 
ferences seen in the MiniBooNE data; the variations of 
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the signal cover the bias introduced through the cho- 
sen unfolding method. The uncertainty on the observ- 
able CC7r + cross section is constrained by measurements 
within the data [iH [46| assuming no bin-to-bin correla- 
tions. The total cross-section uncertainty on the observ- 
able CC7T cross-section measurement is 5.8%. 

The creation and propagation of optical photons in 
the mineral oil is referred to as the optical model. Sev- 
eral ex-situ measurements were performed on the mineral 
oil to accurately describe elastic Rayleigh and inelastic 
Raman scattering, along with the fluorescence compo- 
nents [HI iHjJ. Additionally, reflections and PMT effi- 
ciencies are included in the model, which is defined by a 
total of 35 correlated parameters. These parameters are 
varied, in a correlated manner, over a set of data-sized 
MC samples. The uncertainty calculated from these MC 
samples contains an additional amount of statistical er- 
ror; however, in this analysis, the bulk of the additional 
statistical error is smoothed out by forcing each MC to 
have the same underlying true distributions. The opti- 
cal model uncertainty on the total cross-section measure- 
ment is 2.8%. 

Variations in the detector electronics are estimated as 
PMT effects. The first measures the PMT response by 
adjusting the discriminator threshold from 0.1 photoelec- 
trons (PE) to 0.2 PE. The second measures the correla- 
tion between the charge and time of the PMT hits. These 
uncertainties contribute 5.7% and 1.1% to the total cross 
section respectively. 

The dominant uncertainty on the cross-section mea- 
surements comes from the uncertainty of 7r + — > tt° and 
7r + — > X in the mineral oil occurring external to the ini- 
tial target nucleus where the effects of ir + ir° and 
7r + — > X internal the target nucleus are included in these 
measurements. The uncertainty on the ir + — > tt° cross 
section on mineral oil is 50% and for tt + — ^ it is 35%. 
The uncertainties come from external data [361 ] . These 
uncertainties affect this measurement to a large degree 
because of the much larger observable CC7r + interaction 
rate (by a factor of 5.2). While this uncertainty is small in 
the observable CC-7r + cross-section measurements, many 
CC7r + that undergo either ir + — > tt° or tt + — > X in the 
mineral oil end up in the CC7T candidate sample. The 
uncertainty on the CC7r° cross sections of tt + — > ir° and 
7r + — ► X from observable CC7r + interactions is 12.9%. 
The uncertainty applied to CC(NC)multi-7r, NC7r + , and 
other backgrounds is included in the background cross- 
section uncertainties. 



C. Discussion 

The total systematic uncertainty, from all sources, on 
the observable CC7r° total cross-section measurement is 
18.7%. The total uncertainty is found by summing all 
of the individual error matrices. The largest uncertainty, 
7r + — > 7T° and 7r + — > X in the mineral oil, is 12.9%; the 
flux uncertainties are 10.5%; the remaining detector and 



neutrino cross-section uncertainties are 8.6%. The total 
statistical uncertainty is 3.3%. Table HTT1 summarizes the 
effects of all sources of systematic uncertainty. Clearly, 
the limiting factor on the measurement is the understand- 
ing of 7r + — > 7T° and 7r + — > X in mineral oil external to 
the target nucleus. The two simplest ways to reduce this 
uncertainty in future experiments are to improve the un- 
derstanding of pion scattering in a medium, or to use a 
fine-grained detector that can observe the ir + before the 
charge-exchange or absorption. Beyond that, gains can 
always be made from an improved understanding of the 
incoming neutrino flux. 

TABLE III. Summary of systematic uncertainties. 

source uncertainty 

7T+ ->• 7T° & 7T+ -»■ \ in mineral oil 12.9% 

v„ flux 10.5% 
v cross section 5.8% 
detector electronics 5.8% 
optical model 2.8% 

total 18.7% 



VIII. RESULTS 

This report presents measurements of the observable 
CC7T cross section as a function of neutrino energy, and 
flux-averaged differential cross sections in Q 2 , E^, cos 9^, 
p^o I , and cos 9^o . These measurements provide the most 
complete information about this interaction on a nuclear 
target (CH2) at these energies (0.5-2.0 GeV) to date. 
Great care has been taken to measure cross sections with 
minimal neutrino-interaction model dependence. First, 
the definition of an observable CC7r° interaction limits 
the dependence of these measurements on the internal 
7T + — > 7T° and 7r — > X models. Second, most of the mea- 
surements are presented in terms of the observed final- 
state particle kinematics further reducing the dependence 
of the measurements on the FSI model. Any exceptions 
to these are noted with the measurements. 

The first result, a measurement of the total observable 
CC7T cross section is shown in Fig. El This measurement 
is performed by integrating Eqn. [TJ over neutrino energy. 
This result does contain some dependence on the initial 
neutrino-interaction model as the unfolding extracts back 
to the initial neutrino energy; however, this dependence 
is not expected to be large. The total measured sys- 
tematic uncertainty is 18.7%, and is slightly higher than 
the uncertainties presented for the cross sections with- 
out this dependence. The total cross section is higher 
at all energies than is expected from the combination of 
the initial interaction [301 ] and FSI as implemented in 
Nuance. An enhancement is also observed in other re- 
cent charged-current cross-section measurements [44|,[45j]; 
however, the enhancement is a factor of 1.56±0.26 larger 
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FIG. 8. (color online) The total observable CC7r° cross sec- 
tion as a function of neutrino energy. The uncertainty is dom- 
inated by 7r + charge-exchange and absorption in the mineral 
oil external the target nucleus. The total systematic uncer- 
tainty on the cross section is 18.7%. The central-value mea- 
surement, uncertainties, and correlation matrix are tabulated 
in Table IVII of the appendix. 



than the prediction here. 

The differential cross sections provide additional in- 
sight into the effect of final-state interactions. By neces- 
sity, these cross sections are presented as flux-averaged 
results due to the low statistics of these measurements. 
Except for the measurement in Q 2 , these measure- 
ments are mostly independent of the underlying neutrino- 
interaction model (though there is a slight influence from 
the unfolding). The flux-averaged cross section, differen- 
tial in Q 2 (Fig. is dependent on the initial neutrino- 
interaction model because it requires knowledge of the 
initial neutrino kinematics. The measured Q 2 is un- 
folded to the Q 2 calculated from the initial neutrino and 
final-state muon kinematics. This measurement shows 
an overall enhancement along with a low-Q 2 suppres- 
sion (relative to the normalization difference) with a to- 
tal systematic uncertainty of 16.6%. A similar disagree- 
ment is also observed in the CCn + cross-section measure- 
ment [45| . 

The kinematics of the \x~ are fully specified by its ki- 
netic energy and angle with respect to the incident neu- 
trino beam as the beam is unpolarized. Fig. [TU] shows 
the flux-averaged differential cross section in /i~ kinetic 
energy. Like the total cross section, this shows primarily 
an effect of an overall enhancement of the cross section 
as the fi~ is not expected to be subject to final-state 
effects. The total systematic uncertainty for this mea- 
surement is 15.8%. The flux-averaged differential cross 
section in n~-v^ angle (Fig. [TTj) . shows a suppression of 
the cross section at forward angles, characteristic of the 



FIG. 9. (color online) The flux-averaged (0.5 < E v < 
2.0 GeV) differential cross section in Q 2 with total systematic 
uncertainty of 16.6%. The central-value measurement, uncer- 
tainties, and correlation matrix are tabulated in Table IVTT1 of 
the appendix. 
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FIG. 10. (color online) The flux-averaged (0.5 < E v < 
2.0 GeV) differential cross section in with total systematic 
uncertainty of 15.8%. The central-value measurement, uncer- 
tainties, and correlation matrix are tabulated in Table IVTTT1 
of the appendix. 



low-Q 2 suppression. As the shapes of the data and the 
Rein-Sehgal model as implemented in Nuance are fun- 
damentally different in this variable, the unfolding proce- 
dure required a reduction of the number of bins in order 
to be stable. The total systematic uncertainty is 17.4%. 
The 7T° kinematics yield insight into the final-state in- 
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FIG. 11. (color online) The flux-averaged (0.5 < E v < 
2.0 GeV) differential cross section in cos # M with total sys- 
tematic uncertainty of 17.4%. The central-value measure- 
ment, uncertainties, and correlation matrix are tabulated in 
Table HXl of the appendix. 



teraction effects and are also fully specified by two mea- 
surements: the pion momentum and angle with respect 
to the neutrino beam direction. Fig. [12] shows the flux- 




FIG. 12. (color online) The flux-averaged (0.5 < E v < 
2.0 GeV) differential cross section in |p T n| with total system- 
atic uncertainty of 15.9%. The central- value measurement, 
uncertainties, and correlation matrix are tabulated in Table IXl 
of the appendix. 

averaged differential cross section in |p w o|. The cross 
section is enhanced at low momentum and in the peak, 



but agrees with the prediction at higher momentum. A 
similar disagreement is also observed in the NC7r° cross- 
section measurements [l|. Interactions of both the nu- 
cleon resonance and pions with the nuclear medium can 
cause the ejected tt° to have lower momentum. The total 
systematic uncertainty is 15.9%. Fig. [13] shows the flux- 
averaged differential cross section in ifi-v^ angle. The 




FIG. 13. (color online) The flux-averaged (0.5 < E v < 
2.0 GeV) differential cross section in cos#„.o with total sys- 
tematic uncertainty of 16.3%. The central-value measure- 
ment, uncertainties, and correlation matrix are tabulated in 
Table IXII of the appendix. 

cross section is more forward than the prediction. The 
total systematic error is 16.3%. 

Each of the cross-section measurements also provide 
a measurement of the flux-averaged total cross section. 
From these total cross sections, all of the observable 
CC7T cross-section measurements can be compared. Ta- 
ble IIVI shows the flux-averaged total cross sections cal- 

TABLE IV. Summary of the flux-averaged total cross sec- 
tions calculated from each cross-section measurement. The 
average cross section is calculated assuming 100% correlated 
systematics. The flux-averaged neutrino energy is {E v ) 9 = 
0.965 GeV. 



measurement (er 


, [xl0~ 39 cm 2 ] 


<t{E v ) 


9.05±1.44 


da/dQ 2 


9.28±1.55 


da/OE,, 


9.20±1.47 


da/d cos 6fi 


9.10±1.50 


da /d\p n o\ 


9.03±1.54 


da/d cos 6 w o 


9.54±1.55 




9.20±1.51 



culated from each measurement. The measurements all 
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agree within 6%, well within the uncertainty. While 
all measurements use the same data, small differences 
can result from biases due to the efficiencies and un- 
foldings. The results are combined in a simple aver- 
age, assuming 100% correlated uncertainties, to yield 
(a) 9 = (9.20 ± 0.3 stat . ± 1.5l syst .) x lO" 39 cm 2 /CH 2 at 
flux-averaged neutrino energy of (E v )^ = 0.965 GeV. 
The averaged flux-averaged total cross-section measure- 
ment is found to be a factor of 1.58 ± 0.05 sta t. ± 0.26^4. 
higher than the Nuance prediction. 



IX. CONCLUSION 

The measurements presented here provide the most 
complete understanding of CC7r° interactions at energies 
below 2 GeV to date. They are the first on a nuclear 
target (CH 2 ), at these energies, and provide differential 
cross-section measurements in terms of the final state, 
non-nuclear, particle kinematics. The development of 
a novel 3-Cerenkov ring fitter has facilitated the recon- 
struction of both the ir° and \i~ in a CC7r° interaction. 
This reconstruction allows for the measurement of the full 
kinematics of the event providing for the measurement of 
six cross sections: the total cross section as a function of 
neutrino energy, and flux-averaged differential cross sec- 
tions in Q 2 , Ep, cos Op, \p„o\, and cos 9^o. These cross 

sections show an enhancement over the initial interaction 

i — i. 

model 30] and FSI effects as implemented in Nuance. 
The flux-averaged total cross section is measured to be 
(d) 9 = (9.20 ± 0.3 stat . ± \.h\ syst ) x 10~ 39 cm 2 /CH 2 at 
mean neutrino energy of (E v )^ = 0.965 GeV. These mea- 
surements should prove useful for understanding incoher- 
ent pion production on nuclear targets. 



trino and muon 4-momenta. Fig. [14] shows the 
background-subtracted reconstructed invariant mass for 
data and the MC expectation. The data has not been cor- 
rected for cut efficiencies. The fact that the data peaks 
somewhat below the A(1232) resonance while Nuance 
peaks at the resonance implies the model is not properly 
taking into account final state interaction effects; how- 
ever, it is observed that the MiniBooNE data is almost 
completed dominated by the A(1232) resonance. This 
shift can also be interpreted as an effective change in the 
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FIG. 14. (color online) Invariant mass of the hadronic system 
for both data (points with error bars) and the MC prediction 
(solid line) for the signal mode. The dotted line indicates the 
location of the A(1232) resonance. 
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Appendix: Hadronic invariant mass 

The background-subtracted reconstructed nucleon res- 
onance mass is calculable from the reconstructed neu- 



recoil mass, W, of the hadronic system. Additionally, it 
has been verified that CCQE interactions, which do not 
involve a resonance, peak at threshold (not displayed). 



Appendix: Tables 

The tables presented in this appendix are provided to 
quantify the flux, Fig. [7J and the cross-section measure- 
ments, Figs. E M EH HU ECU and [TJ 
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TABLE V. The predicted flux. Tabulated are the central value, total systematic uncertainty, and the correlation matrix. The bin boundaries are tabulated as the 
low edge of the bin with the last bin giving both the low and high edge. 
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TABLE VI. The total observable CC7r° cross section as a function of neutrino energy (Fig. [8j . Tabulated are the central value (CV), the total systematic error, the 
statistical error, and the correlation matrix for the systematic error. The correlations for the statistical errors are small and not tabulated. The CV, total systematic 
error, and statistical errors are multiplied by 10 39 . The bin boundaries are tabulated as the low edge of the bin with the last bin giving both the low and high edge. 
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TABLE VII. The flux-averaged observable CC-7r° differential cross section in Q 2 over the flux range E v £ (0.5 — 2.0 GeV) (Fig. [9]). Tabulated are the central value 
(CV), the total systematic error, the statistical error, and the correlation matrix for the systematic error. The correlations for the statistical errors are small and not 
tabulated. The CV, total systematic error, and statistical errors are multiplied by 10 39 . The bin boundaries are tabulated as the low edge of the bin with the last bin 
giving both the low and high edge. 
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TABLE VIII. The flux-averaged observable CC7r° differential cross section in over the flux range E v £ (0.5 — 2.0 GeV) in bins of E^ — m M (Fig. [TOjl . Tabulated are 
the central value (CV), the total systematic error, the statistical error, and the correlation matrix for the systematic error. The correlations for the statistical errors 
are small and not tabulated. The CV, total systematic error, and statistical errors are multiplied by 10 39 . The bin boundaries are tabulated as the low edge of the bin 
with the last bin giving both the low and high edge. 
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TABLE IX. The flux-averaged observable CC7r° differential cross section in cos# M over the flux range E v £ (0.5 — 2.0 GeV) (Fig. Tabulated are the central value 
(CV), the total systematic error, the statistical error, and the correlation matrix for the systematic error. The correlations for the statistical errors are small and not 
tabulated. The CV, total systematic error, and statistical errors are multiplied by 10 39 . The bin boundaries are tabulated as the low edge of the bin with the last bin 
giving both the low and high edge. 
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TABLE X. The flux-averaged observable CC7r° differential cross section in |p T o| over the flux range E v £ (0.5 — 2.0 GeV) (Fig. [T^J. Tabulated are the central value 
(CV), the total systematic error, the statistical error, and the correlation matrix for the systematic error. The correlations for the statistical errors are small and not 
tabulated. The CV, total systematic error, and statistical errors are multiplied by 10 39 . The bin boundaries are tabulated as the low edge of the bin with the last bin 
giving both the low and high edge. 

Bin edge [GeV/c] 0.00 0.10 0.15 0.20 0.25 0.30 0.40 0.50 0.60 0.80 1.00-1.40 



CV [xlO 39 cm 2 /GeV/c] 
Stat. 
Total Syst. 



0.00 
0.10 
0.15 
0.20 
0.25 
0.30 
0.40 
0.50 
0.60 
0.80 
1.00-1.40 



4.92 26.65 32.90 28.99 19.02 13.65 7.41 4.27 1.90 0.87 0.19 
±0.37 ±1.10 ±1.05 ±0.87 ±0.61 ±0.37 ±0.31 ±0.33 ±0.25 ±0.35 ±0.21 
±1.34 ±4.94 ±5.00 ±4.31 ±3.09 ±2.49 ±2.01 ±1.14 ±0.63 ±0.40 ±0.79 
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TABLE XI. The flux-averaged observable CC7r° differential cross section in cos^o over the flux range E v £ (0.5 — 2.0 GeV) fFig. 113ft - Tabulated are the central value 
(CV), the total systematic error, the statistical error, and the correlation matrix for the systematic error. The correlations for the statistical errors are small and not 
tabulated. The CV, total systematic error, and statistical errors are multiplied by 10 39 . The bin boundaries are tabulated as the low edge of the bin with the last bin 
giving both the low and high edge. 
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